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NOMENCLATURE 


Note: Many Intermediate angles and locations are defined as they 

are introduced in the text in the section where the deflection 
calculations procedure is described. They are not listed in 
this nomenclature list. 


Letters 

AN 


d 


-2 

e 


d 

P 

D 

a 

f 

Af 


g 

G 



L . 
1 


I, 

m 



P 


P 

m 

P 

o 

r 


r 

a 


Index of refraction of the transparent walls of an 
SFLT tube 

Diameter of LDV beam at location where its intensity 
is 1/e^ of the maximum intensity 

Fringe spacing in the LDV viewing volume 

Diameter of an LDV seed particle 

Diameter of the light collecting aperture 

Focal length of the LDV transmitting lens 

Band width (frequency response) of photodetector 

The constant of gravitational acceleration 

Light scattering gain parameter 

Light scattering gain parameter for back scatter 

Light scattering gain parameter for forward scattering 

Lean limit, in mole percent, of the ith component of 
a fuel mixture 

Lean limit, in mole percent, of a fuel mixture 
Mole percent of the ith component of a fuel mixture 
Atmospheric pressure 
Measured pressure at a rotometer 

Power of the laser beam being used in an LDV system 

Radius location of an undeflected LDV measuring volume 
in an SFLT 

Distance from measuring volume to light collector 
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t 

V 

t 

V 

p 

V 

8 


V 

V 

V 

^wheel 

Ax 


o 
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Greek letters 

6 


8 ' 


n 

q 

0 


X 

W 



T 


Inner radius of the SFLT 
Outer radius of the SFLT 
Tine 

The velocity of the gas that an LDV seed particle Is 
Immersed In 

The velocity of an LDV seed particle 

The final settling velocity of an LDV particle in the 
earth's gravitational field 

Doppler signal visibility, i.e. amplitude of peaks to 
amplitude of pedestal 

Flow rate calculated from a rotometer setting 
Flow rate indicated by a rotometer setting 

Rotational velocity of a calibration wheel 

Displacement of the LDV measuring volume in the x direction 
due tc passage through the SFLT (negative values mean 
motion toward the source) 

Off axis location of the undeflected LDV measuring volume 
Off axis location of the deflected LDV measuring volume 


Angle between the beams that produce the LDV measuring 
volume 

Calculated angle between the beams that produce the LDV 
measuring volume after deflection by the SFLT tube walls 

Quantum efficiency of a photo detector 

Net rotation of the LDV measuring volume due to the 
deflection of the intersecting beams by passage through 
the walls of the SFLT 

Wave length of the laser light used in an LDV 

Viscosity of the gas that an LDV seed particle is 
immersed in 

Doppler frequency for a test of the system 
Density of the gas in an LDV experiment 
Density of an LDV seed particle 

Characteristic particle relaxation time for an LDV seed 
particle to respond to a step change in the gas velocity 
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Subscripts 

o Initial state 

.015 15 mV laser power 

2.0 2W laser power 

Abbreviations and 
Acronyms 

BDFD Beam Deflection Flame Detector 

LDV Laser Doppler Velocimeter or Laser Doppler Velocimetry 

NASA National Aeronautics and Space Administration 

SFLT Standard Flammability Limit Tube 

SNR Signal to Noise Ratio 

STP Standard Temperature and Pressure 

Thermal Systems Incorporated 


TSI 
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ABSTRACT 

The results of the study reported here have bearing on the 
use of the Laser Doppler Veloclmeter, (LDV), to measure the flow 
associated with the passage of a flame through a Standard Flammabil- 
ity Limit Tube (SFLT). There are four major results of this study. 
These are: 

1. It has been shown that using standard ray tracing calcula- 
tions one can predict the displacement of the LDV volume and the 
fringe rotation within the experimental error of measurement. 

2. The flow velocity vector field associated with passage 
of an upward propagating flame in an SFLT has been determined. 

3. It has been determined that the use of a light interruption 
technique to track particles is not feasible. 

4. It has been shown that a 25 mW laser would be adequate 
for LDV measurements in the Shuttle or Space Lab. 
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I 

INTRODUCTION AND SYNOPSIS OF LDV THEORY 


1.1 Introduction 

Gabeous fuel-air mixtures exhibit a range of composition over which a 
flame will propagate spontaneously for long distances from a strong ignition 
source*' These limits are apparatus dependent and they are strongly 
influenced by the earth’s gravity field. Approximately 50 years ago the 

United States Bureau of Mines developed a "standard" flammability apparatus 

which consists of a vertical 51 mn interior diameter tube of approximately 1.8 
meters length. For limit determination mixtures of fuel and air ere placed in 
the tube. One end of the tube is then opened and ignition is effected at the 
open end. If the flame propagates the entire length of the tube the mixture 

is said to be flammable. If the flame extinguishes somewhere in the tube the 

mixture is said to be nonflammable and if the flame does not leave the region 
of the ignition source one is observing an ignition limit, not a flammability 
limit, due to the fact that the ignition energy is not sufficient to ignite a 
flame. It has been observed that in thi.» flammability tube the limit is 
different for upward vs. downward propagation of the flame. In particular for 
the methane-air system the lean limit for upward propagation is 5.25Z CH^ 
while the lean limit for downward propagation 5.86X CH^. Also the shape of 
the downward propagating flame and the mechanism of extinction are markedly 
different than those for the upward propagating flame. 

Because of the large effect of gravity on the nature of flame propagation 
and on the actual flammability limit in a standard flammability tube NASA 
Lewis Research Center has identified this experiment as one which should be 
studied under zero g conditions using either the Shuttle or Space Lab. 
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Preliminary studies to validate and define the usefulness of Shuttle or Space 
Lab experiments have been continuing for the past 4-5 years. Prior to the 
work reported In this report Investigators at the University of Illinois have 
used a single bean Laser Doppler Veloclaetcr systea to aeasure the vertical 
component of the flow velocity vector associated with the passage of the flame 
along the center line of the tube and have used holographic techniques to 
determine the density field associated with the flame and have photographed 
the extinction of the flame using an Image Intensified . All these 
experiments were performed under one g conditions. Additionally the behavior 
of the flame under zero g conditions was studied at NASA Levis Research Center 
using the two second drop towed. In these experiments the shape of the flame 
was photographed and it6 propagation speed through the tube was measured. 

Lean limits could not be determined however because under zero g conditions 
near the lean limit the flame had a propagation speed of only 5 cm/sec and 
thus traveled approximately 10 cm during the 2 second drop time. This 
distance of propagation was Insufficient to say with certainty that 
extinguishment was really a flammability extinguishment and not an Ignition 
extinguishment. 

The purpose of the work that Is reported herein Is to 1) determine the 
effect of the presence of the plexiglass tube on the location of the Laser 
Doppler Veloclmeter intersection volume and on the orientation of the fringes 
in that intersection volume, 2) simultaneously aeasure both the radial and 
the vertical components of the velocity vector during flame passage, 3) 
investigate the feasibility of using an lnterupted particle track technique to 
follow the flow associated with flame passage and 4) assess the feasibility 
of using a 25 nW Laser Doppler Velodteter in the shuttle or Space Lab to 
measure the flow velocity field associated with the passage of a flame under 


zero g conditions. 
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1 .2 The Laser Doppler Veloclaeter 

A Laser Doppler Veloclaeter Is a tool which allows one to sake nonintru- 
slve aeasureaents of the velocity of particles suspended In a fluid such as 
water or air. If these particles are saall enough (dlaneter < 1 micrometer) 
they will follow the notion of the flow quite faithfully even under condltlona 
where relatively high frequency turbulent fluctuations are occurring. The 
principle of the Laser Doppler Veloclneter Is based on the fact that an opti- 
cal laser produces a single frequency beam of coherent light of very high 
intensity and very small diameter. The word coherent here means that the 
light is oscillating In phase with itself just as the 60 cycle line current or 
the nlcrowave signals that are used to transmit information. The principle of 
operation of the Laser Doppler Velocimeter Is quite simple. A single beam of 
coherent light Is split Into two beams and these beams are aligned to be 
parallel to each other and some distance apart. Then these two beams are 
Intercepted by a single lens whose center line lies exactly on the centerline 
of the two beams. This lens focuses the two beams at the same point in space 
producing an intersection volume In which the coherent nature of the light 
causes interference fringes to appear as shown in Fig. 1. 



Figure 1. Representation of fringe pattern created by the 
intersection of two coherent light beams. 
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These stationary fringes ere equally spaced and fora planes that are perpen- 
dicular to the plane of the two Incident beams . Spacing between the fringes 
la controlled by the wavelength of the light and their Intersection angle, 

dj • X sin (S/2) 

where df is the spacing, \ is the wavelength of the light, and S Is the angle 
between the two incident beams. Since the beams themselves have essentially 
Gaussian intensity distribution, the fringe region is actually a small double 
conical region In which the fringes have a Gaussian Intensity distribution. 

Vhen a flow that Is seeded with small particles passes through the beam 
each Individual particle is illuminated by the interference pattern. A sensl 
tlve phcfomultlplier focused on the beam Intersection volume will record the 
Doppler burst shown in Fig. 2. After the pedestal is removed the burst 
appears as shown In Figure 3. Here the overall amplitude envelope is 
essentially Gaussian in time and the frequency , f, of the internal signal 
yields the velocity component of the flow in the direction normal to the 
fringe orientation 



Figure 2- Doppler burst before pedestal removal. 
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Figure 3. Doppler burst after pedestal removal. 

It Is obvious that with this type of arrangement one cannot discriminate 
between forward or reverse flow. This problem is handled by frequency shift- 
ing one of the beams using an optical technique. The frequency of that beam 
is usually shifted about 0.01 to 40 megahertz. This causes the fringes in the 
Doppler volume to travel through the volume at an extremely high velocity, 
much higher than any gas velocity that one would wish to observe. Thus, with 
frequency shifting, the frequency of the Doppler signal is increased mark- 
edly. Electronic instrumentation is used to subtract this frequency off the 
burst before it is analyzed. Using this technique one can discriminate be- 
tween forward and backward flows as they pass through the Doppler volume. 

another innovation two laser beams of different colors are split into 
two pairs of parallel beams whose planes are orthogonal to each other and the 
lens now forms two sets of orthogonal fringes in the Doppler volume. Optical 
processing of the signal allows one to observe each of the two colored fringe 
patterns with two different photomultipliers and thus one can simultaneously 
obtain two components of the flow velocity at the same point in space. 

The work presented in this report was performed using a two component 
frequency shifted Laser Doppler Velocimeter powered by a 2 watt Argon ion 


laser. 
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II 

REFRACTIVE EFFECTS 

2.1 Introduction 

It Is anticipated that either a one color or a two color laser Doppler 
system will be used to measure the flow velocities associated with flame 
passage in a standard flammability tube under zero g conditions. It is 
planned to orient the laser relative to the tube such that the beams will lie 
in vertical and horizontal planes (along the tube and across the tube) and the 
probe volumes will traverse the major axis of the tube that is orthogonal to 
the axis of the two beams. Thus, both the axial and radial components of the 
flow velocity will be measured along that axis. This will require a separate 
experiment for each location of the probe volumes. Unfortunately because the 
tube is circular and has a finite thickness and index of refraction the laser 
beams entering the tube will be deflected and will intersect at a location 
slightly removed from their intended intersection location. Also because the 
two beams enter the tube at different angles relative to the surface of the 
tube when the intersection volume is not at the center line of the tube the 
viewing volume may be slightly rotated because of passage through the plexi- 
glass tube. For a vertical tube and the vertical set of intersecting beams 
the rotation of the viewing volume will not cause a rotation of the fringes 
because the fringes are in a horizontal plane. However, for the horizontal 
beam, the fringes arc vertical and therefore rotation of the viewing volume 
will also cause a rotation of the orientation of the fringes. This means that 
the direction of the horizontal velocity vector that is measured will not be 


strictly radial. 
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2*2 Deflection Celculetlons 

Ray tracing techniques can be used to calculate the effect of the plexl- 
glass tube on all of these properties of the Intersection volume. The calcu- 
lation is reasonably complex and will be outlined step by step. First con- 
sider the case of horizontal beams In a vertical tube. Figure 4 shows the 
upper beam (In the drawing) at an angle alpha to the center line of the tube, 
Intersecting the outer tube wall, of radius R2, at a point A and placed so as 
to Intercept the major axis of the circular tube at point B If no deflection 
were to occur. 



Figure 4. Location of the upperbeam Intersection 
with outer tube wall when undeflected 
beam Is displaced a distance y (point B 
from the center of the tube, 0). 

This point B is a distance y Q above the center of the tube, 0. The angle 
beta, which is the angle between the ray and a normal to the surface of the 
tube where It enters the tube, can be determined by using the sine law: 
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'_o m Z 

•in 0 " »in (a + 90°) 

or 

y 

0 - arcsln (-=^ ain (a + 90°)] 

*2 

Then the triangle ACO in Figure 4 can be used to calculate the location of the 
point at which the beam intercepts the outside wall of the tube. 

( y 2 ) u " ®2 sin + P) 

(*2^ u * *2 COS + P) 


Since the index of refraction of the tube material is higher than that of air 
the entering beam will be deflected in a direction towards the local normal to 
the surface of the tube. In other words the angle delta in Figure 5 will be 
less than the angle beta. 



Figure 5. . A schematic drawing of the index 

of refraction law AN = sin 0/sin 6. 
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The index of refraction law yields the relationship 

6 » arcsln sin fj] 


where AN is the index of refraction. We must now find the location of the 
beam when it intercepts the inner surface of the tube and the angle that it 
makes relative to the normal to the inner surface of tube at the point where 
it exits the tube wall. Figure 6 shows the construction for that calculation. 



Figv-e 6. The location of the intersection of 

the light beam with the inner surface 
of the tube (point D). 


The sine law can be used to determine y' which is somewhat less than 180°. If 
we use the principle value of gamma, y *» 180 - y', since sin y * sin (180 - 
y') we obtain the relationship 
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and we can calculate the angle epsilon using the relationship 

e » y - 6 

We note from Figure 6 that we can now determine the (x, y) location of the 
point at which the beam exits the tube as well as the exit angle relative to 
the y - 0 axis. 

(y^) u * Rj sin (a + $ + e) 

(x. ) * R. cos (a + B + e) 
v 1 'u 1 


We must now calculate the angle to the x axis that the beam makes when it 
exits the tube wall. Referring to Figure 7 the index of refraction law states 
that 

n * arcsin [AN sin y] 

inner 
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and the geometry of Figure 7 shows that the desired angle £ u is given by the 
relationship 

For the lower beam the geometry Is only slightly different. Notice from 
Figure 8 that for the lower beam the angle that was 90 + o for the upper beam 
becomes 9G - a* 



Figure 8. The calculation of the location of the 
point that the beam intercepts the wall 
(point A) for the lower beam. 


However, sin (90 + a) * sin (90 - a). Thus 8 and 6 are equal for the upper 
and lower beam for the same displacement y Q . Inis *eans the y and e are also 
equal for the upper and lower beam. However the location of the intersection 
of the lower beam with the outer wall of the tube must be calculated using the 
equations 
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(y 2 ) L " r 2 ,ln " a) 

(x 2 ) L - 1*2 cos ($ - a) 

its exiting location must be calculated using the equations 

( y l “ R 1 Bln ~ a + £ ) 

(xj) L ■ Rj cos ($ - u + e) 

and its exit angle using the equation. 

C L «g-a+c-n 

Figure 9 shows how the intersection point, the intersection half angle, and 
the rotation of the viewing volume is calculated. 



Figure 9. The geometry used to calculate the actual location 
of the intersection volume. 
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Since C L la negative we can define two ratios 


and 


A 


tan (c u ) 





B 


fr1 


We therefore have two equations in two unknowns 


and 


A fx. ) - Ax, 

v 1 'u c 



y 


o 


b(x 1 ) l **0 * y 0 ~ ( y A 


which can be solved for x Q and y 0 ' to yield 


( y l^L ^ y l ^u + A ^*l ^u + B ^*l ^1 


(A + B) 


and 


y o' m B x ol 


Therefore the displacement of the beams intersection point is given by 


Ax * -x_ 
o 

and 


A y “ y 0 ’ - y 0 
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The half angle between the two beaas la given by 

8 /2 * ('» ■ s .) /2 

and the net rotation by the equation 

9 ■ + «J /2 

A computer program was written to perform these calculations. It is listed in 
Appendix A. A tabulation of these calculated values is given in Table 1 for a 
beam half angle of 10.8° and in Table II for a half angle of 5.515°. Their 
import will be discussed after the calculation for the vertical beams and the 
experimental measurement of deflections are presented. 

TABLE I 

Effect of plexiglass tube on location and rotation 
of the horizontal beam intersection volume 
(for half angle of 10.8°) 

(refer to Figure 9) 


1 

o 

y'o* 00 

Ax, ran 

0, deg 

672, deg 

0.0 

0.000 

0.000 

0.000 

10.80 

1.0 

1.000 

-.001 

-.082 

10.80 

2.0 

2.000 

-.006 

-.165 

10.80 

3.0 

3.000 

-.013 

-.250 

10.80 

4.0 

4.000 

-.024 

-.337 

10.80 

5.0 

5.000 

-.038 

-.426 

10.80 

6.0 

6.000 

-.055 

-.520 

10.80 

7.0 

7.000 

-.077 

-.619 

10.80 

8.0 

7.999 

-.102 

-.724 

10.80 

9.0 

8.999 

-.133 

-.837 

10.80 

10.0 

9.999 

-.170 

-.959 

10.80 

11.0 

10.998 

-.213 

-1.093 

10.80 

12.0 

11.997 

-.263 

-1.240 

10.80 

13.0 

12.996 

-.323 

-1.405 

10.80 

14.0 

13.995 

-.394 

-1.590 

10.80 

15.0 

14.993 

-.478 

-1.801 

10.80 

16.0 

15.990 

-.579 

-2.045 

10.80 

17.0 

16.987 

-.701 

-2.331 

10.80 

18.0 

17.982 

-.851 

-2.672 

10.80 

19.0 

18.975 

-1.037 

-3.088 

10.80 

20.0 

19.964 

-1.275 

-3.606 

10.80 

21.0 

20.948 

-1.586 

-4.273 

10.80 

22.0 

21.921 

-2.009 

-5.170 

10.79 

23.0 

22.876 

-2.622 

-6.457 

10.79 

24.0 

23.786 

-3.604 

-8.520 

10.79 

25.0 

24.564 

-5.631 

-12.840 

10.79 
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y 0 > 1 

0.0 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9,0 

10.0 

11.0 

12.0 

13.0 

14.0 

15.0 

16.0 

i:< .0 

18.0 

19.0 

20.0 

21.0 

22.0 

23.0 

24.0 

25.0 


TABLE II 

Effect of plexiglass tube on location and rotation 
of the horizontal beam intersection volutae 
(for half angle of 5.515°) 

(refer to Pigure 9) 


y'o* ■» 

Ax, an 

0, deg 

872, deg 

0.000 

0.000 

0.000 

5.515 

1.000 

-.001 

-.083 

5.515 

2.000 

-.006 

-.168 

5.515 

3.000 

-.013 

-.253 

5.515 

4.000 

-.024 

-.341 

5.515 

5.000 

-.038 

-.433 

5.515 

6.000 

-.056 

-.528 

5.515 

7.000 

--078 

-.629 

5.515 

7.999 

-.K>4 

-.736 

5.515 

8.999 

-.136 

-.851 

5.515 

9.999 

-.173 

-.977 

5.515 

10.998 

-.217 

-1.114 

5.515 

11.997 

-.269 

-1.265 

5.515 

12.996 

-.330 

-1.435 

5.515 

13.995 

-.403 

-1.627 

5.515 

14.993 

-.490 

-1.847 

5.515 

15.990 

-.595 

-2.102 

5.515 

16.986 

-.722 

-2.403 

5.515 

17.981 

-.880 

-2.764 

5.515 

18.973 

-1.077 

-3.208 

5.515 

19.961 

-1.331 

-3.766 

5.515 

20.942 

-1.668 

-4.496 

5.515 

21.912 

-2.136 

-5.496 

5.515 

22.857 

-2.831 

-6.977 

5.515 

23.742 

-4.014 

-9.498 

5.515 

24.404 

-7.143 

-16.374 

5.515 
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Determining the effect of off-exls behavior on the deflection of the 
vertical beams la a considerably more complex problem than determining the 
deflection of the horizontal beams. In this case since both beams enter the 
tube at the same off-axis location and the same angle to the horizontal they 
will both be deflected In an identical manner (as mirror Images). As with the 
horizontal beams the elevation of the intersection volume (l.e. the z direc- 
tion) will not be displaced. Figure 10 Is a three dimensional drawing Illus- 
trating the geometry used to calculate the refraction of the beam when It 
enters the tube wall. The angle alpha in this figure lies in a horizontal 
plane and the angle 10.8° is the half angle of the two beams and lies in the 
vertical plane. The angle beta lies in a plane which contains the beam and 
the normal to the wall surface at point 0. This plane is tilted by an angle 
gamma to the horizontal plane. This plane also contains the refracted beam 
inside the wall and the refracted beam inside the tube itself. The first step 
in the calculation is to determine the location of the point at which the beam 
intersects the tube wall relative to the center of the tube at z * 0. Refer- 
ring to Figure 11, for a horizontal displacement of the beam y Q the angle 
alpha and the location (x 2 * y 2 > z 2^ can ** calculated using the equations 


a - arcsin (y Q /R 2 ) 
X 2 ■ R 2 cosa 

*2 - y Q 

z 2 “ *2 tan 


1 
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Figure 10 The geometry when a beam in the vertical plane 
enters the tube at an off axis position 


18 



Figure 11. Location of the Intersection point of the beam 
In the vertical plane with the outside surface 
of the tube and the angle, a, that the beam makes 
with that surface in the horizontal plane 


A 
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Now referring beck to Figure 10, If we as suae thet the line eegwent AO has 
unit length we can calculate the angles beta and gams using the following 
foraulas 


L 

L 

L 

L 

L 

Y 

B 


j - tan 10.8 
2 * tan a 



- ere tan (Lj/l^) 


■ arccos 


2 2 
L 4 + L S~ 

2L 4 L 5 


3 


The Index and refraction law and the angle beta may now be used to calculate 
the angle delta which lies In the plane OCED Inside the tube wall. 


6 - arcsln [i— sin fi] 


As stated earlier, this plane Is the same plane as the plane AOB that contains 
beta. However, on this side of the point 0 this plane Is below the horizontal 
plane OFGD. In order to continue the calculation the angle delta nust be re- 
resolved back Into the angles eta and zeta which lie 1:< the horizontal and 
vertical planes, respectively. Pigure 12 Is a drawing of the OCED plane and 
it Illustrates how the value of the angle epsilon and the length Lj ’ are 
calculated using the following equations 


4 

t»> - arcsln [— slnB] 

L 3 

t - ( 180 - « - w) 
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Figure 32. Geometry for the calculation of e and the 

length in the COD plane. Note that in general 
there are no right angles in these triangles. 


Specifically, t 1 '.^ angle omega is calculated using the sin law. Epsilon is 
determined for the triangle OED and is calculated using the sin law in the 
triangle OED. Next notice that the triangle DFC and the included triangle DGE 
shown on Figures 10 and 13 are similar triangles. This allows the calculation 
of Lj* • 

L2 * * l^L-j f /L^ 



C 

Figure 13. The geometry used to calculate • 
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In Figure 14 the triangle OPCD is drawn and used to calculate the value of 
eta* 


D 



Figure 14. Geometry for the calculation of n in the OFGD plane* 

Note that the angle OFD is a right angle in this plane. 


or 


(a - n) ■ arctan (1* 2 - l^'J 


n * a - arctan (L 2 - l^') 


Note that eta Is the angle the refracted beam makes to the vertical plane that 
runs through the normal to the surface- The length of the lines OG and GE are 
then calculated to allow the determination of the angle zeta 


« 1.0/cos (a - n) 


Referring back to Figures 10 and 13 the length Lj * can now be calculated and 
this allows the determination of zeta, which is the angle the refracted b»am 
makes to the horizontal plane. 


X 
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Lj' - L 2 ' tan Y 
C » arctan (Lj'/L^) 


We next need to determine the location where the beam exits the wall to 
enter the tube. Figure 15 shows the geometry associated with this problem. 



Figure 15. Calculations of the angular and (xj , y^) 
position when the beam exits the wall. 


The sin law can be used to calculate the value of the angle theta and simple 
addition can be used to determine the angle alpha* at which the beam exits the 
wall. 


therefore 


and 


R 1 _ *2 w *2 

sin n sin (180 - 0) sin 9 

0 * arcsin t ( Rj /Rj ) sin n] 


a* * a + 9 - n 
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Figure 15 also Illustrates the geometry used to calculate the locations (xj 
and y^) where the beam exits the wall. 

Xj - Rj cos o' 
yj ■ Rj sin a' 

The value of z^ can be determined by calculating the path length of the 
beam in the plexiglass tube wall and noting that its angle to the horizontal 
plane passing through a normal to the surface in the tube is given by the 
angle zeta. 

p * /[(*2 ~ * i ) 2 + ■ yi > 2 ^ 

Zj * z 2 - P tan C 

To determine the exit angle for the beam in the horizontal and vertical planes 
we can use the same procedure as was used to determine the angles eta and 
zeta, only now the entrance angle, 10.8, becomes the angle zeta and the off- 
axis angle normal to the surface, alpha, becomes the angle theta. The geo- 
metry is identical to the geometry used to calculate the behavior of the beam 
when it . ntered the tube wall, except that nes values of eta prime and zeta 
prime are calculated after using the index refraction formula for the exiting 
beam 


6 * 


arcsin (AN sin6') 
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.Figure 16 The geometry that defines the direction of the beam 
as It exits the tube wall and enters the tube. Note 
that the angle 10.8 In figure 9 becomes £ and the angle 
a In figure 9 becomes 0. Also note that 6 > 0. 
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Figure 16 is a three dimensional drawing similar to Figure 10 except that 
It shows that when the beam exits the tube wall to enter the tube the re- 
fracted beam OE exits at an angle to the normal to the surface which Is larger 
than the incident beam AO. This Is because the Index of refraction of the gas 
In the tube Is less than that of the wall. Figure 17 shows the geometry that 
is used to calculate the x location of the beam Intersection point. 



Figure 17. Geometry used to calcuate Ax for the 
vertical beam intersection volume. 

The location of the intersection point x' and the Ax displacement of this 
location may be calculated using the following formula. 

Ax * Xj - x' - Xj - Zj/tan 


i 
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The geometry for calculating the lateral displacement and rotation of the beam 
la ahown In Figure 18. The equatlona for thl8 calculation are aa follows. 

For the beam rotation 

6 - o' - n* 

and for the lateral displacement 

yj m y\ + x i ’ tan ® 

A separate computer program was written to perform the calculations for the 
vertical beams. It is listed in Appendix B. Tables III and IV contain the 
results of these calculations for the vertical beams. Their import will be 
discussed in section 2.4 after the presentation of the measured deflections. 






Figure 18. The geometry used to calculate the angle of 

rotation of the beam (a' - 6') and the actual 
(y') location of the intersection volume. 
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TABLE III 

Effect of plexiglass tube on the location and rotation 
of the vertical beam intersection volume 
(for a half angle of 10.8°) 


*** y'o* ““ 


0.0 

0.000 

1.0 

1.018 

2.0 

2.036 

3.0 

3.053 

4.0 

4.071 

5.0 

5.090 

6.0 

6.108 

7.0 

7.126 

8.0 

8.145 

9.0 

9.164 

10.0 

10.184 

11.0 

11.204 

12.0 

12.224 

13.0 

13.246 

14.0 

14.268 

15.0 

15.292 

16.0 

16.319 

17.0 

17.347 

18.0 

18.380 

19.0 

19.418 

20.0 

20.466 

21 .0 

21.527 

22.0 

22.616 

23.0 

23.759 

24.0 

25.054 

25.0 

28.017 


(refer to Figure 18) 


Ax, ran 


-1.065 
-1.065 
-1.066 
-1.068 
-1.070 
-1.073 
-1.076 
-1.080 
-1.085 
-1.090 
-1.096 
- 1.102 
-1.109 
-1.116 
-1.124 
-1.131 
-1.139 
-1.147 
-1 .154 
-1.161 
-1.167 
-1.170 
-1.170 
-1.165 
-1.152 
-1.126 


0, deg 

B 72, deg 

0.000 

10.80 

-.096 

10.80 

-.194 

10.80 

-.292 

10.80 

-.394 

10.80 

-.499 

10.80 

-.609 

10.80 

-.726 

10.80 

-.849 

10.80 

-.982 

10.80 

-1.126 

10.80 

-1.284 

10.80 

-1.458 

10.80 

-1.652 

10.80 

-1.872 

10.80 

-2.124 

10.80 

-2.416 

10.80 

-2.759 

10.80 

-3.172 

10.80 

-3.677 

10.80 

-4.312 

10.80 

-5.142 

10.80 

-6.281 

10.80 

-7.973 

10.80 

-10.899 

10.80 

-22.760 

10.80 
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TABLE IV 

Effect of plexiglass tube on the location and rotation 
of the vertical beam intersection volume 
(for a half angle of 5.515°) 

(refer to Figure 18) 


o * “> 

y' 0 » ■» 

Ax, mm 

0, deg 

872, deg 

0.0 

0.000 

-1.050 

0.000 

5.515 

1.0 

1.018 

-1.050 

-.095 

5.515 

2.0 

2.036 

-1.051 

-.191 

5.515 

3.0 

3.053 

-1.052 

-.288 

5.515 

4.0 

4.071 

-1.055 

-.388 

5.515 

5.0 

5.089 

-1.057 

-.492 

5.515 

6.0 

6.107 

-1.061 

-.601 

5.515 

7.0 

7.126 

-1.065 

-.715 

5.515 

8.0 

8.145 

-1.070 

-.837 

5.515 

9.0 

9.163 

-1.075 

-.968 

5.515 

10.0 

10.183 

-1.081 

-1.110 

5.515 

11.0 

11.203 

-1.087 

-1.266 

5.515 

12.0 

12.223 

-1.094 

-1.438 

5.515 

13.0 

13.245 

-1.102 

-1.631 

5.515 

14.0 

14.267 

-1.109 

-1.849 

5.515 

15.0 

15.291 

-1.117 

-2.098 

5.515 

16.0 

16.317 

-1.126 

-2.387 

5.515 

17.0 

17.346 

-1.134 

-2.727 

5.515 

18.0 

18.378 

-1.142 

-3.136 

5.515 

19.0 

19.416 

-1.149 

-3.637 

5.515 

20.0 

20.463 

-1.155 

-4.269 

5.515 

21.0 

21.524 

-1.159 

-5.094 

5.515 

22.0 

22.611 

-1.160 

-6.228 

5.515 

23.0 

23.753 

-1.156 

-7.914 

5.515 

24.0 

25.045 

-1.144 

-10.832 

5.515 

25.0 

27.998 

-1.119 

-22.685 

5.515 


2.3 Deflection Measurement 

To check the calculations, experimentF were performed to photograph the 
beams as they were refracted by the tube. A 40 cm section of tubing with 
properties identical to those of the Standard Flammability Limit Tube (SFLT) 
was used for the experiments. It was a transparent polymethylmethacrylate 
tube, having an Inside diameter of 50.72 nn>, a thickness of 2.79 mm, and a 
— — Jtmf ractlva tiff 1«*9» 
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Because Che refractive effecta were avail in relation to the correa- 
ponding amount of off-axis positioning, accurate alignnent of the tube aection 
with the laser beams was necessary. Figure 19 depicts the arrangenent for 
refraction measurements of the horizontal beams. A 79 cm x 137 cm table which 
translated in two dimensions was used to act as a moving support for the tube 
section and photographic system. The laser beam source and transmitting len6 
were fixed so that -vovement of the table in the Y-direction v >ld move the 
tube for off-axis me. ements. 

Both the horizontal and vertical beams were used for alignment. The 
table surface wa6 first carefully leveled to form a horizontal plane parallel 
to and about 3 cm below the plane formed by the two horizontal beams. The two 
axes of table translation were aligned along the X and Y axes as defined in 
Fig. 19. Alignment was accomplished without the SFLT in place by positioning 
a mirror mounted 90 degrees to the table surface at the beam crossing, facing 
the transmitting lens as shown in Figure 20. When the table was properly 
oriented, the reflected beams coincided with the Incident beams and 
translation of the table along the Y axis did not alter beam coincidence. 
Translation of the table in the negative X direction (toward the lens) caused 
the four reflected beams to form a perfectly centered rectangular pattern 
inside the location of the Incident beams on the lens surface. 

A 90 degree draftsman triangle was used to check the table alignment with 
respect to the vertical beams, as 6hown in Fig. 21. When the table wa6 
properly aligned, the vertical edge of the triangle would ju6t "touch" both 
beams simultaneously as the table was translated to the right. 

Once the table was aligned, a flat front surface mirror was laid face up 
on it directly below the beam crossing. The 40 cm long SFLT section was 
mounted directly above the mirror. The tube was adjusted to be perfectly 
vertical by aligning it with its reflection in the horizontal mirror. 
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Figure 20 Alignment test of laser beams with table. 

Drawing depicts table tilted slightly forward. 



figure 21 Arrangement for table alignment with vertical 
beams. 
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The vertical beams were blanked off and the horleontal beans were 
photographed with a Canon A-l 35 ■& canera with a Vivltar 210 ■» telephoto 
lens and #1 close-up lens attachment. The canera was nounted on a tripod 
approximately 10 cn above the tube top, as shown In Fig. 19. To position the 
camera to look directly down the tube axis, the reflection of the camera lens 
In the horizontal mirror was centered in the camera field of view. This 
oriented the film plane In the camera normal to the tube axis. The camera was 
then moved horizontally until the reflection of the lens lay exactly in the 
center of the image of the tube bottom. This procedure centered the camera 
directly over the tube axis. 

The circular inside edge of the tube bottom served as a reference in 
making measurements from the photographs. Since the horizontal beams were 
only 1 cm above the tube bottom, the camera system operating at f8 provided 
adequate depth of field. The distance from the camera film plane to the beams 
was 50 cm, and to the tube bottom, 51 cm. Parallax effects in the measure- 
ments were accounted for. 

To check flatness of field, a fine grid was photographed with the above 
photographic system from a distance of 50 cm, and the negative was measured on 
the Eausch & Lomb lCx optical comparator. The grid spacing was found to be 
constant in the 50 mm diameter region of concern. 

The alignment procedure for the measurement of horizontal refraction has 
been described. The procedure for vertical refraction measurements was simi- 
lar because the vertical position y* and deflection angle 8 could be 
measured using the set up shown in Fig. 19. The only difference was that the 
camera had to be mounted horizontally to view the tube from the -Y side to 
determine values of Ax. To avoid obtaining a refracted image of the already 
refracted beams, a 60 degree arc centered on the -Y axis was cut out of the 
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tube and measurements were aide of the beams on the +Y aide of the tube* 

To conduct the measurements, the beans nuat be centered in the tube with 
high precision. The transnlttlng optics were first translated In the y 
direction of Fig. 19 until all secondary bean reflections formed a symmetric 
pattern on the tube walla- This located the probe volume on the x axis to an 
accuracy of ± 0.06 mm Iti the y direction. The table was then translated in 
the positive and negative x directions until the refracted and reflected 
horizontal converging beams formed two single spots on the opposite inside 
walls- The table position was measured with aetrlc dial gages , and the tube 
center was assumed to be located exactly at this point. Accuracy of centering 
the probe volume on the x axis with respect to the y axis was estimated at t 
0-5 mm. A water vapor spray at the base of the tube was used to make the 
beams visible for photography. Double exposures were made: one of the beams 

in the tube center for reference and another at one of the selected off-axis 
positions • 

The negatives of the photographs were examined on a Bausch & Lorab 10X 
optical comparator- Accuracy in the measurements of beam crossing position 
are estimated to be t 0.4 sd, and of the probe volume rotation, ± 0.4 
degrees. The accuracy was better at smaller off-axis positions. 

The results of the experiment are tabulated with the calculated values at 
a number off-axis positions in Table V,for the case when Pj/2 * 5. 15°. Only 
this half angle was used for the measurements because an examination of the 
calculated results, tables I through IV, showed that the refraction effects 
were largeet for this half angle. Off-axis distances do not have integer 
aetrlc values because the translating table was calibrated in Inches, and 
because small corrections in tran/iated position were made using the reference 



TABLE V 


Comparison of Measured and 


Calculated Refraction Results 


* 1/2 


5.515“ 


A. Horizontal beams 


Off axis 

y ' 

7 o 


Ax 


6 


Position 

mm 


ran 


degrees 


■m 

calc* 

mess 

calc* 

meas 

calc* 

meas 

5.50 

5.50 

5.4 

-0.05 

0.0 

-0.48 

-0.5 

10.89 

10.89 

10.7 

-0.21 

-0.2 

-1.10 

-1.2 

16.22 

16.21 

16.0 

-0.62 

-0.6 

-2.17 

-2.2 

19.32 

19.29 

19.3 

-1.16 

-1.2 

-3.39 

-3.3 

21. 5A 

21.45 

21 .3 

-1.92 

-1.9 

-5.04 

-4.6 

22.51 

22.39 

22.3 

-2.49 

-2.4 

-6.25 

-6.3 

23.65 

23.44 

23.2 

-3.60 

-3.7 

-8.62 

-8.5 

25.30 

- 

24.0 

- 

-5.8 

- 

•12.5 

Vertical beams 

(y’o and 

e ) 





Off axis 

y’o 




9 


Position 

nun 




degrees 


mm 

calc** 

meas 



calc** 

meas 


5.50 

5.60 

5.4 

-0.55 

-0.6 

10.89 

11.11 

10.7 

-1.25 

-1.1 

16.22 

16.54 

16.2 

-2.46 

-2.2 

19.32 

19.75 

19.4 

-3.84 

-3.5 

21.54 

22.11 

21 .7 

-5.67 

-4.8 

22.51 

23.19 

22.8 

-7.09 

-5.7 

23.65 

24.59 

24.2 

-9.79 

-8.2 

ertical beams 

(Ax only) 





Ax 


Position 

mm 


mm 

calc** 

meas 

0.00 

0 

m 

p-H 

1 

-1.1 

5.36 

-1.06 

-1.1 

10.71 

-1.08 

-1.1 

19.29 

-1.15 

-1.2 

21.43 

-1.16 

-1.2 

22.51 

-1.16 

-1 .0 

23.57 

-1.15 

-1.1 

24.64 

-1.13 

-1.0 


Linearly interpolated from table II 
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beams In the photographs. Because the caacia was Mounted In a different 
position to Measure the x-value for the vertical bean Intersection point, a 
separate set of off axis values were used and the off-axis positions were Dot 
exactly equal to those of the other tests > 

2.4 Comparison of Predicted and Measured Refractive Effects 

In all the measurements, as well as in the calculations, the intersection 
angle was essentially unchanged for both sets of beams. The majority of 
measurements agree with the predictions to within the accuracy of 
measurement. Small systematic discrepancies, probably due to errors in tube 
positioning, or to lnhomogenelties in the tube, become pronounced only at 
lt-rge off-axis positions. Overall, out to 22.5 mm, the greatest discrepancies 
between calculations and measurements are only 1.4 degrees in rotation and 0.4 
am in position. The agreement for the horizontal beams is much better showing 
maximum differences of only 0.4 degree and 0.24 am. 

Beyond 22.5 mm, the larger disparities are possibly due to positioning 
errors, since the calculations show that in this region the probe volume 
position and rotation are very sensitive to position changes. Taking into 
account the work described above and the resulting agreement between analytic 
and experimental data, it can be assuoied that the calculations can be used to 
accurately predict the position of the LDV beam intersection volumes. T .t 
should be noted that the rotation of the vertical beams through the angle 
0 does affect the fringe orientation by that anwunt while the rotation of the 
horizontal beams has no effect on the fringe orientation because the 
interference planes will remain horizontal. 
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III. 


VELOCITY MEASUREMENTS USING THE LDV 


3.1 Introduction 

The purpose of this portion of the study Is to simultaneously measure the 
radial and vertical velocity components of the transient flow velocity at 
points located on the tube axis during flame passage for upward propagation of 
a near-lean-limit flame in a standard flammability tube using a two color 
frequency shifted laser Doppler velocimeter. The apparatus, technique and 
results will be presented in some detail in this section. 

Before describing the individual components of the experiment in detail, 
a brief overview of a typical experiment is in order. A standard flammability 
tube (SFLT) meeting the specifications of Coward and Jones * was filled from 
the top with a fuel/air mixture. The fuel and air flows were monitored by two 
separate rotometers and mixed just before entering the tube. The fuel went 
directly from the rotometer to the mixing point, whereas the air, after 
exiting the rotometer, went to a particle generator to collect seeding 
particles needed for LDV measurements. The tube was purged with a volume of 
mixture equal to ten times the volume of the tube. This ensures that the 
metered mixture conposition is the actual composition in the tube. After 
filling, the SFLT was isolated and then opened at the bottom and ignited at 
the bottom using a hot-wire ignitor. 

The gas velocities associated with the upward propagating flame were 
measured using a two watt, two-color, frequency-shifted LDV system. Frequency 
shifting permitted measurement of the flow reversal that occurred when the 
flame passed the point of measurement. Photodetectors in the system converted 
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Che Doppler burets to electrical signals of the saoe frequency, and these were 
fed Into s Counter-type signal processor which converted each electrical 
signal to a voltage proportional to the signal frequency. The analog voltage 
output was recorded on a Tektronix Type 533A four channel oscilloscope. A 
referencing system was constructed to coordinate the time-dependent voltage 
with the changing position of the flame. 

3.2 Flammability Apparatus 

The flammability tube constructed for this study was copied from Levy^ 
and is depicted s c hem." ti tally in Figure 22. It was constructed of 50.80 ms (2 
in) ID, 3.175 mm (.125 in) thick clear plexiglass tubing, and was 1.8 m 
(70.9 in) in ] jiig n, thus meeting the requirements specified by Coward and 
Jones. * 

At each end of the tube a round plexiglass flange 24 cm in diameter was 
affixed eccentrically. The center of each flange served as the center of 
rotation for a plexiglass disk. The upper rotating disk (see Fig. 23) 
contained a plenum chamber where the air and methane used for filling the tube 
were introduced and mixed. There was also an open hole in the disk which 
could be positioned over the tube for venting and cleaning between tests. 

When neither of these two positions were used, the top end of the SFlT was in 
the closed position. 

The rotating disk on the lower flange had a 7 mm diameter hole which was 
positioned under the tube to allow the gas mixture to flow out during the 
filling process (see Figure 24). After filling, the lower disk was rotated to 
the closed position for several minutes to allow the mixture time to become 
quiescent. It was then rotated to the open position for ignition. A heated 
nichrome wire was used as an Ignitor. 
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3.3 The Combustible Mixture 

The fuel used for this study was natural gas. This source was aore 
readily available and much aore inexpensive than bottled gas. An analysis of 
the gas Is presented In Table VI. It was chosen because it Is prlaarlly 
ae thane and ae thane Is the fuel of choice for the proposed space lab 
studies. Since the objective of this experlaent was prlaarlly the assessaent 
of an experlaental technique rather than the interpretation of the actual 
velocity aeasureaents or the exact determination of a limit, there was no need 
to use pure methane. 

The addition of small amounts of inerts to a flammable mixture has been 
shorn by Jones and Kennedy^ to have the same effect on the lean limit as the 
addition of air. Their measured values of flammability limits of mixtures of 
fuel were in excellent agreement with those predicted by Le 'Chateller 's 
rule. Using values of the lean flammability limits from Coward and Jones*, we 
will apply Le 'Chateller 's rule to obtain the lean limit of the mixture. Thus, 

n 

mixture lean limit, L^ « 100/ ( E p. /L. ) 

1-1 

» 100/(88.75/5.3 4 4.65/3.0 + 0.94/2.2 + 0.16/1.9 
4 0.06/1.2 4 0.04/1.5) 

- 5.30 2 

The air and fuel were regulated using Matheson rotoaeters, model numbers 
605 and 603, respectively. They were calibrated by measuring the tine 
necessary at fixed rotameter settings to displace a known volume of water. 

The errors in calibration were estimated to be less than 32. Errors in 
reading the rotometer settings were estimated at less than 22, so an error no 
more than a 52 was expected. 
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The lean Halt wee first determined without the particle dispenser In the 
feed lines and found to be 5.04Z. This la significantly lower than the value 
of 5.30Z determined from Le Chateller's rule. The analysis presented In Table 
VI , however, was 15 months old at the time of testing, so an increase In the 
percentage of ethane or propane might account for the lower limit observed. 
Nvertheless, the lean limit could be measured to an accuracy of at least 
±.05Z. 


TABLE VI 

Natural Gas Analysis* 




Lean 


Volume 

Flammability 

Constituent 

Percentage 

Limit** 

Methane 

88.75 

5.3 

Nitrogen 

4.88 

Inert 

Ethane 

4.65 

3.0 

Propane 

0.94 

2.2 

Carbon Dioxide 

0.50 

Inert 

Butane 

0.16 

1.9 

Hexane-Plus 

0.06 

1.2 

Pentane 

0.04 

1.5 

Helium 

0.02 

Inert 

Hydrogen 

tr 

— 

Oxygen 

tr 

— 

Argon 

tr 



100.00 



* Analysis from Illinois Power Company 

** Values from Coward and Jones ( 1) 

3.4 Particle Seeding for the LDV 

A laser doppler velocimeter measures velocities of particles which are 
seeded in the flow. The requirements in this experiment were that the 
particles must: 1) accurately follow the flow; 2) scatter enough light to be 
detected; 3) be numerous enough to provide an adequate data rate; 4) keep 
their integrity in the hostile flame environment; and 5) absorb a minimum 
amount of heat. It is therefore desirable to use small but highly reflective 
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particles at a concentration as high aa possible without significantly 
affecting the flaae temperature. 

Particle response to changes in fluid velocities is characterised by the 
particle relaxation time: 

t - o d^/ 18 u 
P P 

The particle velocity after a step change in fluid velocity (at a later time 
t) is then given by 

, v -t/t 

... - (» - V >. 

1 o r o 

The particle settling velocity due to a balance between gravitat ioc«»i forces 

o 

(l/6t d^ g(p p - Pj)) and Stokes drag (3it d p uv g ) is 

v g " d p g(p p ~ p f )/18y 

Reuss found aluminum oxide (AI 2 O 3 ) particles are the best particles 
to use in a flammability tube. They are unaffected by the flame, have a 
specific gravity of 3.8, are roughly spherical, and are readily available from 
commercial suppliers in a variety of sizes from 0.1 microns to 40 microns. 

0.1 micron diameter particles were used in Reuss 's LDV system. Because 
aluminum oxide tends to cake it is impossible to generate a monodlspersed 
powder of this diameter. The agglomerates formed from the particle generator 
were assumed to be less than 4 microns in diameter, 60 the relaxation time in 
the mixture at STP would be no more than 185 microseconds, and the settling 
velocity less than 1.8 mm/ sec. 

A particle generator constructed by Reuss, Figure 25, was used in this 
experiment with one alteration: approximately 52 of the lower chamber volume 

was filled with 2 mm diameter glass beads. Also a portion of the air flow 
entered the dispenser near Its base and collected particles in the fluidized 
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bed cheaber. The elr flow held the beade in twirling suspension, minimizing 
particle caking on the chamber walla. A 1.25 an hypodermic needle with four 
oppoalng jeta Inpinging at lta termination Halted the aaxlaua agglomerate 
else when the flow left the fluidized bed chaaber. After pasting through a 
nettling chaaber, the particle laden air panned through the outlet at the top 
and wan trenaferred via 6 ■& plastic tubing to the alxlng "T" at the top of 
the SFLT plenum. The particle generator caused a back pressure at the exit of 
the air rotaaeter of 20.7 kPa. The flow rate was adjusted by applying the 
relationship $ ■ v/p /P '• 

B B 

There was no neter to measure the actual particle density. However, two 
observations can give clues as to the relative density. First, the particles 
present heat sinks to lower the flame temperature and thus raise the lean 
limit. The lean limit of the seeded mixture was aeasured as 5.07Z, .032 
higher than the unseeded limit. This difference Is not entirely reliable, 
however, because It lies with the above mentioned error range of A .052. The 
second clue Is the flame color. An unseeded flame has a dull violet color, 
while a highly seeded flame appears bright orange-yellow. The seeded flames 
observed in these experiments had a mixed violet and weak yellow color. This 
classification is, admittedly, very subjective; but it indicates that the 
flames were not highly seeded. The flame shape was not visibly altered by the 
presence of the particles, thus it was assumed that they were not 
significantly intrusive. 

3 . 5 The Laser Doppler Velocimeter (LDV) 

The SFLT and laser Installation are shown aounted on a table with a 
translating platform in Figure 24. The table that the LDV is aounted on can 
be translated in two orthogonal directions and x, y. The (x, y) position of 
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Supply 


Air 

Rotameter 



Figure 26 Overall arrangement of experimental apparatus. 
Referencing system is not included. 
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th« Manuring volume could be located In the SPLT to vlthln ±0.001 lnchea. 

Figure 27 deplete an LDV arrangement for forward ecatter collection for 
clarity of lliuetratlon. In theae experlaenta backacatter collection was 
used, and the receiving optica minus the receiving lena were mounted between 
the aecond Bragg cell and the tranaalttlng lena. The tranealttlng lena then 
doubled aa a receiving lena. 

The laser was a Spectra-Physics Model 164-06 two watt argon-ion laser. 

It was powered by a Spectra-Physics Model 265 exciter. The exciter provided 
current to the plasaa tube solenoid and controlled the ion discharge In the 
plasma tube so that a constant laser power output was aalntalned. The power 
of the two watt laser is highest at wavelengths of 488 na (35Z of total power) 
and 514.5 na (40Z of total power), corresponding to the colors blue and green, 
respectively. The beam leaving the laser was first collimated so that the 
probe volume fringes were parallel and so the beam "waist*' or minimum 
diameter, would occur at the point of velocity aeasureaent. The bean was 
resolved into its component colors by a TSI Model 901 dispersion prism. Plane 
mirrors directed the green beams down the main horizontal optical axis and the 
blue beams down a parallel axis to bypass the green beam splitter. The beam 
displacer realigned the blue beam along the main vertical optical axir>. 
Polarization rotators were used to rotate the polarity of the beams 
perpendicular to the plane of the beamsplitter. In this way, the beams were 
resolved into two parallel components of equal intensity. The green beams 
were separated by 50 am in the horizontal plane, and the blue beams by 50 am 
in the vertical plane. Bragg cells were used to shift the frequency of one 
beam of each color and thus allow near zero and negative velocity 
measurements. They were powered by two TSI model 985 LDV frequency shifters, 
each with a shifting range of 10 KHz - 40 MHz. Three pairs of achromatic 
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.Figure 27. LDV assembly for forward scatter collection. 

(Note that while backscatter collection was 
used, a diagram of the forward scatter \ i •> 1 H e 
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lease* of focel lengths 120 m, 250 si, end 600 were available to focus the 
laser beoas. 

The constant of proportionality between the maaured Doppler frequency 
and the flow velocity is the fringe spacing, d f . The green beam fringe 
spacing for the 120 me lens was determined by measuring the output frequency 
from a source of known velocity. A rotating disk was positioned so that the 
outer edge had a maximum downward velocity at the beam crossing. Tiny 
osparaties on the disk surface acted as particles to generate Doppler signals 
at a sufficiently high data rate. The disk diameter was measured as 22.02 ±01 
cm, and the rotational speed was measured with a calibrated strobe as 1784 rpm 
± 2 rpm. This implied an outer edge velocity of 20.57 m/sec ± 0.03 m/sec. 

The Doppler frequency was measured as 15.20 MHz, so the fringe spacing was 
simply calculated as df * v wheel / v n * 20.574/15.20 x 10^ “ 1.354 x 10“ 6 m ± 
.003 x 10“^m. This is within 1.5Z o i the manufacturer's specifications. 

Beam colllmatlon can be accomplished by making the fringes parallel. 

This was checked by translating the rotating disk along the long axis of the 
probe volume and adjusting the colllmatlon until the frequency output (and 
thus the fringe spacing) was nearly constant. After final colllmatlon, the 
fringe spacing varied by less than 0.8Z across the piobe volume where the data 
rate was at least 1/e^ (141) of the maximum rate. The length of the probe 
volume as defined by the.l/e^ data rate positions was measured as 0.67 cm. 

This is of the same order of magnitude as the calculated value of 0.32 cm. 

The receiving optics utilized a color splitter with a dichroic mirror to 
pass green light and reflect blue. The receiving modules focused the 
scattered light from the particles in the probe volume onto the 
photomultipliers. Each receiving module had a narrow band-prrr filter; one 
passed only blue light, the other only green. This virtually eliminated noise 
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from other light eourcee such me the flame. TSI Model 962 photomultipliers 
were centered above each of the receiving modules. Clear apertures of 0.127 
mm (.005 in) were centered between the photomultipliers and the receiving 
modules to block extraneous light not scattered from the probe volume. The 
photomultipliers each had a quantum efficiency of 22X, a sensitivity of 40,000 
V/lumen, and a frequency response of 100 MHz. They converted any impinging 
light signal to an electric signal of the same frequency. 

The output from each photomultiplier was directed to the respective 
frequency shifter for downmixing. With the 120 nm lens, the maximum frequency 

« _ . _ _ _ j _ i j _ . . , _ j i _ «.i v * - - A £ . 

umC to poL ti.ci.co miuv xttg uomuwoiu xa me tuoc w da ca^rp l UAitaaLc ajt u*u rmt* ilic 

next highest frequency shift available was 1.0 MHz, so the green probe volume 
fringes were shifted downward at this value to maximize measurement 
sensitivity. Likewise, the maximum expected frequency due to radial 
velocities in the negative (outward) direction was 0.15 MHz, so the radial 
beams were shifted outward by 0.5 MHz. The frequency shifted outputs were 
transmitted to TSI Model 1990 Counter Type signal processors. They had high 
pass filters to eliminate the pedestal frequency, and low pass filters to 
reduce high frequency noise. The processors also rejected bursts which lacked 
a minimum number of cycles (from particles far from the probe volume 
center.) Best results for the measurement of vertical velocities were 
obtained with the filters set at 0.3 MHz and 2.0 MHz, and a minimum of 8 
cycles per Doppler burst to be measured. For radial measurements, the filters 
were set at 0.3 MHz and 1.0 MHz, also with 8 cycles to be counted. 

The counters produced an analog voltage proportional to the measured 
Doppler frequency. The voltage output was held until a new Doppler burst was 
encountered. The output when displayed on an oscilloscope therefore was a 
series of step functions. If the data rate was high enough, the output 
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•pprosched • ■mooth curve. The proportionality constants of the counters. In 
V/MHz, were checked by input t lug sloe waves from a signal generator. The 
output was always within 0.4Z of the expected voltage. 

3 . 6 Referencing Systen 

The output from the counter was a time-varying voltage proportional to 
the velocity at the point of measurement. It was necessary to reference the 
velocity with the changing position of the flame to make any meaning of the 
data. This was accomplished by determining the flame position at one Instant 
of time and measuring the fieuue speed. 

A Beam Deflection Flame Detector (BDFD) was copied from von Lavante^. 

a 

The beam from a Spectra Physics model 124A 15 mWatt HeNe laser was split by a 
prism assembly into, two parallel, horizontal components. The beams were 
directed through the tube center by making sure the beam reflections on the 
tube wall coincided with the beams themselves. 

Two photodetectors (Motorola MRD 510-7727) were positioned to intercept 
the beams 1.7 m on the opposite side of the tube. The photodetectors were 
powered by a 5 volt DC power supply. In the absence of a flame, the laser 
beams struck the photodetectors, the circuits were completed, and a +5 volt 
deflection was registered on the oscilloscope. The density gradient 
associated with the preheat zone of the propagating flame first deflected the 
lower of the two beams and started an Industrial Timer Corp. digital clock 
(.001 sec resolution). The deflection of the upper beam stopped the clock. 
Single 1.2 tm apertures were placed in the beams on either side of the tube to 
increase the sensitivity of the system. The DC coupled output of the 
photodetector circuit was displayed on the oscilloscope to provide a temporal 
reference mark. The flame velocity was simply calculated by dividing the 
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■assured distance between the two beaas by the clock tiae. 

The exact position of the visible flaae with respect to the lower beam at 
the tine the beam was refracted had to be deternlned photographically. It was 
also desirable to photograph the flaae to determine Its shape for use In a 
velocity profile construction. A Canon A-l 35 mm SRL camera was positioned 
facing the tube with the film plane parallel to the vertical plane formed by 
the two reference beams. When the flame deflected the lower beam, a voltage 
pulse was Input into a Tekronlx Type 555 oscilloscope. The oscilloscope was 
set for external triggering on a negative slope, so the voltage drop instantly 
triggered the oscilloscope and a +20 V gate pulse. The pulse opened a silicon 
control rectifier in a solenoid power circuit developed by Meagher ^ which 
activated a solenoid attached to the shutter release and exposed the film. 

There was a mechanical delay from when the +20V gate pulse was initiated 
to the time of shutter opening. To determine this delay time, a calibration 
scheme developed by Reuss was used. The lens of the Canon A-l camera was 
removed and the camera back was opened. The beam from the 15 milliwatt laser 
was spilt by a prism into two parallel beams spaced 21.1 nm apart, and both 
beams were directed through the camera film plane. The beams struck 
photodetectors to complete two 5 Volt d.c. inputs to the oscilloscope. The 
shutter of the Canon A-l consists of two shades which move from right to left 
when the camera is triggered. The first shade exposes the film as it moves, 
and the second shade covers the film. The movement of the first shade 
permitted the beams to pass through the camera body and complete the circuit. 

To conduct the time delay tests, the oscilloscope was triggered manually, 
so that the sweep began at the same time the gate signal was triggered. 

During the time delay, the beam was blocked by the camera shades and the 
oscilloscope registered zero. The oscilloscope registered a step rise of 5 
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volt* When the first shade uncovered the beaas. The oscilloscope traces were 
photographed on Polaroid Type 667 flln and are presented In Figure 27. 

The results of these experiments indicated that the delay time was 
constant regardless of shutter speed. Photographs at every shutter speed from 
1/125 sec to 15 sec showed a delay tine of 73.3 nsec ± 1.3 nsec for the first 
shade to uncover half of the film plane. The shutter speed was measured as 
3.1 m/sec, or 3.1 nm/msec. 

A photograph of the particle-seeded flame is presented in Figure 28. The 
lower beam is the triggering beam, and the upper beam is from the two watt 
laser to provide vertical reference. Because the upper beam does not lie in 
the camera image plane, parallax effects make it appear slanted. However, it 
passes through the tube center, so the apparent beam separation along the 
tubeaxls served as the reference distance. The actual beam separation was 
5.98 cm. 

Five photographs were taken, and the negatives were examined on the 
Bausch & Lomb 10X optical comparator. Precise measurements were made 
difficult because o; : the ambiguity in the location of the flame's leading edge 
and the effect of smearing due to flame movement. During the exposure, the 
flame traveled (23.42 cm/sec) x 1/125 sec = 0.187 cm. The measured thickness 
of what appears to be the reaction zone (the relatively bright leading edge) 
is approximately 0.25 cm, but this may be in error because the upper edge is 
hard to define. However, near the center of the reaction zone, the 
luminescence is maximum and constant in the vertical direction before 
gradually fading in the upward direction. The leading edge of the visible 
flame at the time of shutter opening was assumed to be at the upper end of 
this zone of maximum brilliance. Because of the intensity decrease in either 
direction, it is safe to assume that the flames leading edge has been located 
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Figure 27 Traces of camera shutter delay test. 
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relative to the location of the Bean Deflection Plane Detector to within ± 

0.15 cb. The refractive effects of the tube on the flane Image were accounted 
for by photographing a grid placed in the tube at the caaera focus. The 
results are tabulated In Table VII. 

Three different photodetector aperture arrangements were used to 
determine if the aperture size had an effect on the time of beam deflection. 
The results, given in Table VIII, indicate that the beam is triggered slightly 
earlier when a larger aperture is used. This is surprising because it was 
assumed that the smaller aperture make the system more sensitive. However, 
not enough tests were performed to determine a reliable correlation. The 
combination of a 1.2 mi aperture on the entrance side and an 0.8 us slit on 
the exit side was easiest to use, so it was decided that the visible flame tip 

TABLE VII 

MEASUREMENTS TO DETERMINE 
REFRACTIVE EFFECTS OF TUBE ON FLAME IMAGE 

r (true) r (measured) 

(mm) (mm) 


0.00 

0.00 

3.00 

2.98 

6.00 

6.00 

9.00 

9.03 

12.00 

12.02 

15.00 

15.00 

18.00 

17.87 

21.00 

20.73 

24.00 

23.21 

Tube internal radius ■ 

25.335 mn 

Measurement accuracy ± 

o 

• 

0 
u> 

1 
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TABLE VIII 

FLAME POSITION AT TIME OF BEAM DEFLECTION 
AS A FUNCTION OF APERTURE ARRANGEMENT 


Aperture - Tube 
Entrance Side 

Aperture - Tube 
Exit Side 

Distance of Flame 
Tip Below Beam 

1.2 an 

circular 

1.2 


circular 

.203 

cm 

1.2 ■» 

circular 

1.2 

mn 

circular 

.215 

cm 

1.2 mn 

circular 

0.8 

Bf» 

slit 

.108 

CD 

1.2 BD 

circular 

0.8 

mn 

slit 

.096 

cm 

0.8 ■& 

slit 

0.8 

nm 

slit 

.064 

cm 


should be considered to be 0.10 cm below the lower triggering beam at the time 
of deflection. 

3.6 Interfacing and Testing of System Components 

Successful backscattering LDV measurements require first adapting the 
system to the given test conditions to produce best results. The output 
signal Is optimized when the data rate and signal to noise ratio are 
niximized, and the data scatter (for laminar flow) is minimized. It was found 
that the quality of the signal was strongly dependent on the following: 
transmitting lens used, the laser power, photomultiplier alignment and 
part icle concentration. 

The calculated probe volume dimensions for the three lenses and green 
beam (514.5 nm) are shown below in Table IX. 
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TABLE IX 

Probe volume diameter d f and length 1 B 
for lenses with different focal lengths 


Lens 

Focal Length 

f 

l . 

120 on 

.0603 ml 

• 317 mm 

230 ml 

.131 ML 

1 • 30 mm 

600 mm 

.302 ML 

7.36 mm 


The 600 mm lens yielded a measuring volume much too large for these tests. It 
was originally felt that the 250 mm lens would be better suited than the 120 
un lens because the larger probe volume size would permit a lower seeding 
density to achieve the sane data rate. However, the reverse was found to be 
true. The greatest source of noise is due to diffuse reflections off the tube 
wall caused by the beams striking the tube. With the shorter focal length 
lens, the bright spots on the wall were more widely separated, and the noise 
was reduced. Thi6 allowed the gain to be increased, thereby making a higher 
data rate possible. 

The performances of the Tracker and Counter were compared. The analog 
voltage outputs were of the same quality for the measurement of radial 
velocities, but the Counter slightly outperformed the Tracker for vertical 
velocity measurements. The data rate was higher and the scatter was less. 

This was surprising, because a floating narrow bandpass filter normally makes 
the Tracker better adapted to noisy signals. On the other hand, the Tracker 
was easier to use because the data rate could be directly monitored. This was 
important for the alignment of the photomultipliers. 

Photomultiplier alignment is the most critical factor involved in 
obtaining a good signal. By removing the photomultiplier, a TSI model 10096 
eyepiece was placed over the receiving module for focusing and rough 
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alignment. When scattered reflect Iona fro* the tube were observed to croaa 
over the probe volume, the signal waa later observed to be noisy. Reflective 
interference was serious for both the vertical and horizontal beaas at r < 

5 an and r < 20 ■*. In these ranges, it waa found advantageous, especially 
for radial measurements, to set the receiving aodules slightly out of focus so 
that the reflected Inages lay outside the 0.127 an dlaaeter area of the 
photomultiplier aperture. This compromise resulted in a significantly 
increased data rate and signal to noise ratio. 

Frequency shifting caused the particles that were settling in the tube 
before the combustion tests to generate a Doppler signal. This permitted 
photomultiplier alignment by finding a peak in the data rate and signal to 
noise ratio. The signal to noise ratio for these purposes was defined as the 
normal data rate divided by the data rate obtained when one beam was 
blocked. Between r * 5 an and r * 20 mm, the signal to noise ratio was always 
et least ten. 

An indication of the accessibility of the LDV at different axial 
positions is the relative data rate obtained during the particle settling 
stage. The data rates from the horizontal and vertical beams were measured as 
a function of radial position, and showed a large dropoff when r exceeded 20 
mm. Horizontal measurements at any given radius were more susceptible to 
reflective interference than vertical measurements, and were made difficult or 
Impossible at r < 3 mm and r > 22 mm. Vertical measurements were possible at 
all axial positions, but were more difficult in the same two regions. 

For reasons unexplained, the blue beams yielded a poor signal compared to 
the green beams for velocity measurements. In an effort to obtain the best 
possible performance, the green beams were used to measure both the vertical 
and radial components. The problem was certainly due to the particular LDV 
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system being used, and In no way auggeata Chat simultaneous data could not be 
collected by a different LDV ays tea. 

The only obstacle to alaultaneous data collection la the fact that the 
two probe voluaes are not coincident. Refractive effects of the tube on the 
radial and vertical probe voluaes cause this problea (see Tables I-V). The 
differences In the y-posltlons of the two aeasurlng voluaes la tolerable at 
all axial positions except very near the tube wall. On the other hand, the x- 
poslt Ions always differ from each other by approxlaately 1 am. The horizontal 
beams can be "forced" to Intersect near the vertical beams by passing them 
through a glass plate before striking the tube. The shift in the probe volume 
location Is approximately 1/3 the glass thickness in the minus x direction. 

This method was tested and determined to be viable If the glass plate can 
be oriented parallel to the transmitting lens to within two or three 
degrees. Gross misalignment causes refraction errors which are greater than 
the attempted probe volume shift. 

The beams can quite easily be located in the probe volume center with 
high precision. The transmitting optics were first translated in the y 
direction until all secondary beam reflections formed a symmetric pattern on 
the tube wall. This located the probe volume on the x axis to an accuracy of 
± 0.06 m In the y direction. The table was then translated in the positive 
and negative x directions until the converging beams formed a single spot on 
the front and back inside wails. The table positions were measured with 
metric dial gages, and the ^ center was assumed to be located exactly 
midway between the two readings. Accuracy of centering the probe volume on 
the x axis with respect to the y axis was estimated at ± 0.5 mm. 
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3.7 Experimental Procedure end Keeulte 

It has been observed by Reuss^ and othera that the speed at which a near 
lean-1 ialt flaae propagates up an SFLT Increases with test nuaber if the tube 
is not cleaned between runs. The problem can be eliainated by cleansing the 
tube wall with water after every run. In these tests, the tube walls were 
first wiped with dry chaaols to reaove any particles clinging due to 
electrostatic charges. The tube was then washed with a water filled sponge. 

It is important in the drying procedure to generate as saall a static charge 
as possible on the walls to prevent seed particles from clinging to the vall6 
during the tube filling process. Drying the tube entirely with an air flow is 
much too slow a process. Instead, the tube was partially dried with a cotton 
cloth, and then totally dried with an air jet. The experimental procedure for 
a test cycle is given in Table X. 

3.8 Experimental Results 

Tests were run on a 5.272 mixture flame for vertical velocities at r ■ 0, 
3, 10, 15, 18, 20, 22, and 23.5 mn, and for radial velocities at r “ 5, 10, 

15, 18, 20, and 22 mm. At least two tests were performed at each station to 
determine repeatability. TeBts were al 60 performed on a 5.512 mixture for 
vertical velocities at r • 10 and 18 om>. Only one test was made at each 
location for these latter runs. 

The raw data from the 5.272 tests is presented in Appendix C. It 
represents the data which was most closely repeated in comparison tests, and 
for which the flame was observed to be not significantly asymmetric. On the 
figures in Appendix C, the upper straight line is a zero velocity reference, 
representing the voltage corresponding to the amount of frequency shifting. 

The lower straight but broken line Is the d.c. coupled output of the beam 
deflection flame detector. The voltage jump begins when the flame deflects 
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the lover reference beta, end ends when the upper beta Is deflected. The bees 
eeper et ion vie 11.30 cs. The probe voluse ves located at the ease height as 
the lover bees, ao the uppersoat tip of the flase vaa 0.10 cs below the probe 
voluse at the tlae of thla deflection. 


TABLE X 

EXPERIMENTAL PROCEDURE 


• Fill tube with particle laden mixture for 150 sec 

• Close top and bottom flanges simultaneously 

• Align photomultipliers 

• Adjust gain for reasonable data rate 

• Align beam deflection flame detector 

• Open bottom flange and ignite with hot wire 

• Start oscilloscope trace when flame is 0.1 m below probe 
volume 

• Photograph reference trace for zero velocity 

• Vent tube 

• Remove particles from tube wall with dry chamois 

• Wash tube with water filled 6ponge 

• Partially dry tube with cotton cloth 

• Dry tube completely with air jet 

The traces presented in Appendix C were analyzed with a Hewlett Packard 
model 9864A Digitizer. The two voltage scales of the oscilloscope were 
calibrated, and two third order polynomials were determined to account for 
scale nonlinearity. A program was written to convert the voltage data 
directly to velocities with respect to a flame of average propagation speed. 
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The average flame velocity of the 5.27Z mixture vae 22.89 ca/aec with a 
standard deviation of 0.55 cm/sec, and that of the 5.51Z mixture* 23.80 
ca/aec, with a standard deviation of 0.45 ca/aec. 

The reduced data la presented In Tables XI through XIII. It was obtained 
by first curve-fitting the output data from the digitiser, then aeasurlng 
velocities from the curves at regular Intervals. Vertical velocities for the 
5.27Z aixture at r ■ 0, 5, 10, 18, 20, and 23.5 an are the average of two 
separate tests. z*0 Is defined at the uppermost tip of the flame. The shape 
of the flame was measured on the optical comparator for use In the velocity 
profile construction. A composite of two symmetric flames is drawn in Figure 
29, with velocity vectors drawn in a coordinate system which moves with the 
flame. Refer to the photograph of Figure 28 for comparison to the actual 
flame. 

The accuracy of velocity measurements is largely dependent on the scatter 
of data. The experimental scattter in this investigation was usually no more 
than ± 1.6 cm/sec for vertical velocities and ± 0.9 cm/sec for radial 
velocities. Since turbulence is assumed to be absent, a mean curve is 
considered to represent the true velocity. Ambiguity of the placement of the 
mean curve was approximately one half the maximum data scatter, so the 
measurement accuracy was extimated at ± 0.8 cm/sec for vertical velocities, 
and ± 0.45 cm/sec for radial velocities. 

Vertical velocity measurements are possible at all axial positions, but 
the signal is degraded by reflective interference at r<3 mm and r>22mm. 

Radial velocity measurements are possible in the region 3 mm <r< 22mm. 

Signal dropout was found to occur in radial velocity measurements at r ■ 
18 mm during flow reversal. This may be due to the beams becoming temporarily 
uncrossed or to a particle density gradient near the flame zone. It is not 
unreasonable, however, to interpolate on the range where the di olnts are 
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Table Xi Reduced vertical velocity data for 5.27 % mixture. Velocities are in cm/sec, positive downward, 
relative to upward traveling flame of V p = 22.89 cm/sec. Z is vertical position in tube, 
positive downward, referenced from flame tip. 
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Table XI i Reduced radial velocity data for 5.27% mixture. Velocities are in cm/sec, positive 
radially inward. Z is vertical position in tube, positive downward, referenced 
from flame tip. 
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Table XII i Velocity data for 5.51% mixture. Velocities are in cm/sec. Vertical 
velocity is positive downward, relative to upward traveling flame of 
V p = 23.80 cm/sec. Radial velocity is positive radially inward. Z is 
vertical position in tube, positive downward referenced from flame tip. 




Figure 29. Gas velocities in a coordinate system that is moving 
with the flame. No radial velocity information was 
obtained at r * 23.5 mm. The lengths and direction 
of the arrows represent local velocity and direction 
of the flow at the point where the vector originates. 
The scale is represented by the vector length on 
the left side of the figure. 
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IV 

INTERRUPTED PARTICLE-TRACKING 
(AN ALTERNATE METHOD) 

4.1 Introduction 

It was shown in Section 2 that the refractive effects of the curved tube 
walls have a deleterious effect on the position of the LDV Intersection 
volume. To avoid these problems, an Interrupted particle-tracking technique 
was considered as an alternate means of measurement. Particle-tracking was 
assessed as being feasible because it was known to have some advantages over 
LDV. For Instance, particle-tracking allows velocity measurements at all off- 
axis positions in a single experimental run. 

The basic theory behind interrupted particle-tracking is quite simple. 

As in LDV, the goal is to measure velocities of minute particles which are 
seeded in the flow. The particles are illuminated by a high intensity light 
source which Is interrupted at a fixed frequency for a fixed fraction of 
time. In this way, a moving particle appears on a photograph as a series of 
short 6treaks. By knowing the frequency and interruption fraction and by 
measuring the distances between and length of successive streaks, particle 
velocities may be determined. 

4.2 Experimental Arrangement 

In this study, the two watt laser again served as the light source. To 
measure vertical and radial velocity components, a vertical plane which 
bisects the tube and is perpendicular to the viewing axis of a photographic 
camera was illuminated by a vertical plane of laser light. Such a vertical 
plane of laser light was produced by directing the laser beam through a 
horizontally mounted circular glass rod, as depicted In Fig. 30. The rod 
spreads the light so that a vertical plane of light is formed which passes 




Figure 3 0 


Circular glass rod is parallel to X-axis and spreads 
laser light in a Y-Z plane through the tube main axis. 
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through the tube. To produce a strobe effect, s rotating toothed wheel was 
placed in the path of the laser bean ahead of the glass rod. The illumination 
frequency could be varied by changing the number of teeth on the wheel or Its 
rotational speed. 

A Canon A-l camera was positioned facing the tube so that the film plane 
was parallel to the plane of light illuminating the particles, and was focused 
on this plane. Alignment of the camera film plane and the light plane was 
accomplished in ways similar to those described in section 2. The 
photodetector circuitry used in the LDV study to create a reference time mark 
was used In these tests to trigger the camera automatically. Details of the 
triggering system are given in section 3.5. When the flame deflected the 
lower beam, there was a 73.3 msec delay to the time of shutter opening. 

4.3 Photographic Limitations and Optimization 

A photograph of the particles during flame passage showed many different 
series of streaks, but no flame image 6ince the flame swept across the field 
during the exposure. The flame position at the time of shutter opening was 
determined by measuring the flame speed and referencing the flame with respect 
to the lower triggering beam by applying the shutter delay time. It was 
important to identify the first streak of every particle and assign it a time 
t«0. There is an error that has a maximum value equal to the time when the 
particles are not illuminated, since the particles might not have been 
illuminated at the time of shutter opening. To ensure that most particle 
trails started within the photographed region, the exposure needed to be 1/15 
sec or less. Otherwise, most particles were observed to enter the photograph 
at one edge and leave at another, leaving no clue as to the time of passage. 



70 


Such abort exposures required a aet of photographs to record the entire 
velocity history of the flame passage, and hence defeated one goal, that of 
obtaining all velocity data from a single flame. 

A related problem was encountered by Levy 4 . He attempted to use the 
flame Itself as reference by limiting his exposure to 1/100 sec. However, the 
particles barely moved during such a short exposure, and little Information 
could be obtained. With longer exposures, the flame Image was smeared, and 
referencing was Impossible. 

It was necessary to find a combination of particles and photographic 
system which would yield clear images of particle streaks. The aluminum oxide 
particles used for LDV experiments were highly reflective and were considered 
the best suited for these experiments. To determine the best photographic 
system, te6ts were conducted by filling the tube with particle laden air and 
photographing the particles under strobe illumination. It was immediately 
observed that a portion of the laser light that was used to illuminate the 
particles also remained in the tube walls because of internal reflection. 

This portion of the light that was reflected inside the tube when the light 
sheet first struck the tube, and again when it exited the tube eventually was 
radiated from the tube and obscured particles at off-axis positions greater 
than 60Z of the tube radiu6. 

This problem can be alleviated by passing the light sheet through a 0.8 
did wide 6 lit before it enters the tube. The slit helps to shield the tube 
from scattered laser light from the dispersed beam so that the tube is 
illuminated only be light which is needed for particle illumination. This 
allowed particle to be seen at off-axis positions up to 70Z of the tube 
radiu6. 

In the determination of the optimum photographic system, there were three 
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parameters which had to be determined: the camera lena ayatem, the film, and 

the developer. Since it was desired to use particles as small as possible to 
follow the flow and minimize heat absorption, a fast photographic aystem ia 
required. Four lens systems for close-up work using the Canon A-l were tried, 
and are listed here in order of greatest speed: 1) Canon 50 am fl.8 lens with 

20 am extension tube; 2) Vivltar f2.8 macro len6; 3) Vivitar 210 mm telephoto 
len6 with Vivitar #2 close-up lens; and A) the Canon 50 am fl.8 lens mounted 
in tandem with a 12 mm extension tube and a Vivitar 2x teleconverter. Four 
different film/developer combinations were tested and are listed here in order 
of greatest speed and resolution: 1) Kodak Recording Film 2A75 (ASA 1000), 

Kodak DK-50 developer; 2) Kodak Tri-X pan film (ASA A00), Optimate-2 
developer; 3) Kodak Tri-X Pan film, (ASA A00), Kodak Microdol-X developer, 
and; A) Kodak Plu6-X film (ASA 100), Kodak Microdol-X developer. When using 
the fl.8 lens and 20 mm extension tube with Kodak Recording Film 2A75, 
aluminum oxide particles as small as 3 microns in diameter (assuming no 
agglomeration) were identifiable in photographs of the upward propagating 
flame. 

Although 3 micron diameter particles could be detected, contrast between 
the particles and the internal reflections was not enough to permit reliable 
identification at different time periods. Particles as large as 9 microns 
diameter were necessary for sufficient contrast, although further refinements 
in experimental technique would have certainly lowered this minimum size. 

Since traces of the smaller particles were difficult to distinguish, positive 
identification of the same particle after an interruption of the beam had 
occurred was sometimes dubious. Thl6 made velocity determination using this 
technique a tedious and time-consuming task. 
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The 15 laser waa tea ted agalnat the 2 W laaer aa the light aource to 
order to predict the performance of the Spacelab laaer. The reaulta obtained 
ualng the two lasers were comparable. The particles were not as well 
Illuminated by the smaller laser, but the tube luminescence was also 
diminished, so that the contrast between the particles and tube was about the 
same. 

Finally, scales should be placed on either side of the tube in the same 
plane as the sheet of laser light after flame passage and photographed as a 
double exposure to provide a scale for the measurements. Particle-tracking 
tests were terminated before this step was implemented. 

4.4 Appraisal of Interrupted Particle-Tracking 

The four main advantages of LDV over particle-tracking are: 1) smaller 

particles may be used; 2) velocity data at greater off-axis positions may be 
obtained; 3) data acquisition and velocity determination are less tedious, 
and; 4) most importantly the interpretation of each interrupted trace on a 
photographic record i6 difficult, and it would be particularly difficult to 
determine the exact relationship between particle position and flame motion. 
Thus, particle-tracking as applied to flammability tube work is probably not 
feasible because it is more vulnerable to error and misinterpretation than 
originally anticipated and it should not be investigated further. 


V 

PREDICTED PERFORMANCE OF A 25 oM SYSTEM 
Power limitations on Spacelab require the use of a 15 irirf (nominal) 
laser. A new 15 wU He-NE laser emits 25 mW at A - 632.8 nm. Two component 
measurements must be made, by means of polarization separation. A 15 aW air 
cooled Argon ion lev^r emits 9 uW at A » 514.5 nm, and 15 ok at A ■ 488.0 
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na, so that color separation can be ueed to reeolve the two component e. 

The following analyale Is based on a two color Argon Ion laser. It 
should be noted that a He-NE laser aay be preferred because the entire output 
power Is concentrated at one wavelength. This Is an advantage when asking 
one-coapooent measurements, such as on the tube axis. 

Ref. 8 defines the signal to noise ratio as 


SNR - 4X10 U (^) 

a 

For identical testing conditions, q , d , Af , ard f are equal. D a /r a is 

q p a a 

assumed to be equal to 1 if a 120 tm collecting lens is used. From Adrian and 

Fingerson (1) visibility V is dependent primarily on the ratio d p /d f> so V 

2 -2 

should also be equal. Then, SNR ~ P^Gd e . The vertical velocity 
measurements made in this experiment required a maximum of 0.24W in each beam, 
and only 0.17W per beam in the range 5 mm <r< 22mm. Assuming the 514.5 nm 
beams are to be used to measure vertical velocities, the power required to 
obtain the same SNR achieved in these tests is: 
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No information was found in the literature on the scattering parameters of 

AI 2 O 3 as a function of collection angle. The particles are not truly 

spherical, and Mie scattering theory aay not accurately apply. However, since 

the ratio d p /df I s estimated to be on the order of 3, the ratio Ggs^FS 

q 

should be quite small, and an estimate of 1/100 is considered conservative. 


In that case, 
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- 10.2 mW required in each beam, or 20.4 a M total. A 
similar analysis for the blue beams indicates a requlrment of 25.4 at A - 
488.0 nn. The above analysis is Intended only for an order of magnitude 
estimate, and indicates that a 15 mW laser operating in forward scatter might 
possibly achieve the same SNR obtained with the 2 W system in backscatter. A 
more sophisticated prediction for the expected performance can be computed by 
a program written by Meyers^, at NASA Langley, which is available at NASA 
Lewis. 

It should be mentioned that any calculation of relative performance may 
be significantly inaccurate. For instance, it i6 difficult to predict the 
reduction in reflective noise attained by off-axis positioning of the 
collection optics, or the scattering parameters of A^O^ as a function of 
collection angle. Perhaps the mo6t reliable method is to compare the vertical 
velocity measurements of this experiment to those obtained by Reuse with a 15 
otf system and use the comparison to predict the performance of the smaller 
laser in radial measurements. 

Reuss measured the vertical velocities in a glass tube with forward 
scatter collection to an accuracy of ±0.5 cm/sec, as compared to ± 0.8 cn/sec 
attained in this experiment. There is reason to believe that Reuss* particle 
concentration was higher, since the particle presence raised the lean limit 
from 5.272 to 5.372, a 0.102 increase, versus the 0.032 increase observed 
here. The higher data rate would allow the gain to be reduced, resulting in a 
greater SNR and less scatter. It is therefore assumed that the two systems 
would yield similar results under Identical testing conditions. Reuss* 
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attempts to measure radial velocities mere hindered by signal dropout for much 
of the time. He attributed this to the beams becoming uncrossed. However, 
beam deflection should only occur for a short time near the reaction tone, 
where density gradients are most severe. Moreover, the radial measurements 
made in this experiment did not exhibit prolonged dropout. Reuss' system was 
not equipped with frequency shifting, and the Tracker most probably lost the 
signal when the velocity was below 5 cm/sec. The radial velocities were found 
here to be below this value for most of the time. His data did show a sharp 
velocity peak near t«0 of about 13 cm/sec, and another more gradual peak near 
t«0.2 sec of approximately 10 cm/sec. This is consistent with the results 
obtained here, and indicates that it was the absence of frequency shifting, 
not beam deflection, which prohibited radial measurements. 

Based on the above comparisons, it is concluded that a 15 li frequency 
shifted two component LDV system operating in forward scatter should achieve 
similar performance to the system used here. There is also the possibility of 
superior performance through optimal positioning of the collection optics and 
transducers. 


VI 

RESULTS AND CONCLUSIONS 

As was mentioned in the introduction, the purpose of this study was to: 
1) calculate and verify by experimental measurement the effect of the 
presence of the circular walls of an SFLT on the location and rotation of the 
fringe patterns of an LDV, 2) measure the radial and vertical components of 
the velocity vector during flame passage, 3) investigate the feasibility of 
using an interrupted particle track technique to follow the flow associated 
with flame passage and 4) assess the feasibility of using a 25MW LDV in the 
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Shuttle or Space Lab to determine the flov velocity field aaaoclated with the 
passage of the flame under aero g conditions. 

All of these tasks have been accomplished. The comparison of the 
calculated and measured location and orientation of the LDV measurement volume 
is quite satisfactory and shows that one can determine the deflection caused 
by the circular tube walls using simple ray tracing techniques. The flow 
associated with flame passage for upward propagation of a near limit flame in 
an SFLT in the earth's gravity field has been measured. The use of the 
interrupted particle track technique for Investigating the flow associated 
with the passage of a transient flame in SFLT has been evaluated and found to 
be inferior to the LDV method. Thus, in this report, it i6 recommended that 
this technique should not be considered again. Finally an evaluation of the 
feasibility of using a 25 oN LDV for transient flow velocity measurements in 
either the Shuttle or Space Lab leads to the conclusion that the use of a 25mW 
LDV in this operation would be feasible. 
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APPENDIX A 


PROGRAM ABC (INPUT , OUTPUT) 

ALPH=10. 8*3.14159/180. 

ANIT=3. 14159/2. 

DO 10 1=1,26 
AI=I-1 

BETAU=ASIN (AI*SIN (ANIT+ALPH)/28 . 15) 

BETAL=ASIN (AI*SIN (ANIT-ALPH) /28 . 15) 

Y2U=28 . 15*SIN (BETAU+ALPH) 

X2U=28.15*COS (BETAU+ALPH) 

Y2L=28 . 15*SIN (BETAL-ALPH) 

X2L=28 . 15*COS (BETAL-ALPH) 

DELU=ASIN( (1.0/1. 49) *SIN (BETAU) ) 

DELL=ASIN( (1.0/1. 49) *SIN(BETAL) ) 

GAMU=ASIN (1 . 125*SIN (DELU) ) 

GAML=ASIN (1. 125*SIN (DELL) ) 

EPU=GAMU-DELU 

EPL=GAML-DELL 

Y1U=25.36*SIN (BETAU+EPU+ALPH) 

X1U=25 . 36*COS (BETAU+EPU+ALPH) 

Y1L=25 . 36*SIN (BETAL+EPL-ALPH) 

X1L=25.36*C0S (BETAL+EPL-ALPH) 

ETAU=ASIN (1. 49*SIN (GAMU) ) 

ETAL=ASIN (1 . 49*SIN (GAML) ) 
ZETAU=BETAU+EPU-ETAU+ALPH 
ZETAL=BETAL+EPL-ETAL-ALPH 

X0=( (Y1L-Y1U)+X1U*TAN(ZETAU)+X1L*TAN(-ZETAL) )/ 
* (TAN (ZETAU) +TAN (-ZETAL) ) 

Y0P=Y1L+ (X1L-X0) *TAN (-ZETAL) 

DELX=-X0 
DELY=YOP-AI 
ANG= (ZETAU-ZETAL)/2.0 
ANG=ANG*180./3. 14159 
AROT= ( ZETAU+ ZETAL ) /2 . 0 
AROT=AROT* 180 ./3. 14159 
PRINT 100, A I , Y0P,DELX, AROT,ANG 
100 FORMAT (5X,F5.1,4F10.3) 

110 FORMAT (5X, "ERROR") 

10 CONTINUE 
END 
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PROGRAM ABC (INPUT , OUTPUT) 

CON-180. /3. 14159 
ANIT-3. 14159/2.0 
ALPH-10. 8*3. 14159/180. 

Z2=28.15*TAN(ALPH) 

DEL*0. 0645223 
Z l*Z2-3 . 175*TAN (DEL) 

Y2*0.0 
TH*0 . 0 
ZE=10.8 
Y0P=0 . 0 

DX=25 . 36- (Zl/ (TAN (ALPH) ) ) 

PRINT 180 , Y2 ,Y0P ,DX ,TH , ZE 
180 FORMAT (5X,F5.1,3F10.3,F10.3) 

DO 10 1=1,25 
AI = I 

ALPH=10. 8*3. 14159/180. 

ALPH1-ASIN (AI/28 .15) 

X2=28 . 15*COS (ALPH1) 

Y2=AI 

Z2=X2*TAN (ALPH) 

DO 10 J=1 , 2 

IF (J.EQ.2)ALPH=ZETA 

IF (J.EQ.2)ALPH1=THETA 

PRINT 150 , ALPH , ALPH1 , ZETA , THETA ,Y2 

ELI -TAN (ALPH) 

EL2=TAN (ALPR1) 

EL3=SQRT (EL1**2+EL2** 2) 

EL4=SQRT(1.0+EL1**2) 

EL5=S<JRT ( l . 0+EL2** 2 ) 

BETA=ACOS ( (EL4**2+EL5**2-EL3**2)/(2.0*EL4*EL5) 
IF (EL2.EQ.0.0)GAM=ANIT 
IF(EL2.EQ.0.0) GOTO 5 
GAM=ATAN (EL1/EL2) 

5 CONTINUE 

IF(J.E0.1)DEL»ASIN( (1.0/1. 49)*SIN (BETA) ) 

IF (J .EQ. 2) DELIAS IN (1 . 49*SIN (BETA) } 

AL=ALPHl*CON 

BE=BETA*CON 

GA=GAM*CON 

DE=DEL*CON 

110 FORMAT ( 5 X," ERROR") 

120 FORMAT (4F10. 5) 

PRINT 140, ELI, EL2,EL3,EL4, BETA 
150 FORMAT (5E15. 5) 

ARG= ( (EL4/EL3) *SIN (BETA) ) 

IF (ARG.GT. 1.0) PRINT 130, AT 
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IF(ARG.GT.1.0)OMEG«ANIT 
IF (ARG.GT. 1.0) GOTO 15 

130 FORMAT ( 5X , " ARG GREATER THAN 1.0, AI*",I3) 
OMEG*ASIN (ARG) 

15 CONTINUE 

EPS-ANIT* 2 . 0-DEL-OMEG 
EP*EPS*CON 

EL3P*EL5* (SIN (DEL) /SIN (EPS) ) 
EL2P=EL3P*EL2/EL3 
ETA=ALPH1-ATAN (EL2-EL2P) 

EL6=1 . O/COS (ALPH1-ETA) 

ELIP=EL2P*TAN (GAM) 

ZETA=ATAN (EL1P/EL6 ) 

ET=ETA*CON 

ZE=ZETA*CON 

IF (J ,EQ. 2)GOTO 20 

THETA=ASIN ( 1 . 143*SIN (ETA) ) 

PRINT 160, ZETA, THETA 
160 FORMAT (2E15. 5) 

ALPHP= ALPH 1-ETA+THETA 
ALP=ALPHP*CON 
Xl=25.36*COS (ALPHP) 

Yl=25 . 36*SIN (ALPHP) 

P=SQRT( (X2-X1)**2+(Y2-Y1)**2) 

Zl=Z2~P*TAN (ZETA) 

OM=OMEG*CON 
TH=THETA*CON 
GOTO 10 
20 CONTINUE 

ALAMB-ALPnp-ETA 
ALA=ALAMB*CON 
X1P=Z1/ (TAN (ZETA) ) 

DX=X1-X1P 
DY=X1P*TAN (ALAMB) 

Y0P=YI-DY 

PRINT 150 ,EL1 ,EL2 ,EL3 ,EL4 ,EL5 
*>RINT 150 , X2 , AL ,BE ,GA ,DE 
PRINT 150,OM,EP,EL3P,EL2P,ET 
FRINT 150 ,EL6 ,EL1P , ZE ,TH , ALP 
PRINT 150,X1,Y1,Z1,Y2,Z2 
PRINT 150 ,ALA ,X1P ,DX,DY,Y0P 
PRINT 1 8 0 , Y 2 , Y0P , DX , ALA r ZE 
PRINT 100, AI ,X1,XIP,Y1,Y1P ,ALA , ZE 
140 FORMAT (4X,5F8. 3) 

10 CONTINUE 
100 FORMAT (7F8. 3) 

END 
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ORIGINAL PACE IS 

Op POOR QUALITY 


O.SV/div 


zero 

veloc 1 ty 



Vertical velocity record 


0.2 sec/div 



Radial velocity record 


0.2 sec/div 


Figure C-2. Velocity measurements at r * 5.0 m. 
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Figure C-4. Velocity measurements at r ^ 15.0 mm. 
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